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Abstract. The double pulsar system PSR J0737— 3039A/B is a highly relativistic double neutron star (DNS) binary, with a 
2.4-hour orbital period. The low mass of the second-formed NS, as well the low system eccentricity and proper motion, point 
to a different evolutionary scenario compared to other known DNS systems. We describe analysis of the pulse profile shape 
over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this 
system, PSR J0737— 3039A, to have a low misalignment between its spin and orbital angular momentum axes, with a 68.3% 
upper limit of 6.1 °, assuming emission from both magnetic poles. This tight constraint supports the idea that the supernova 
that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core. 
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INTRODUCTION 

The PSR J0737-3039A/B system H H has provided 
pulsar astronomers and neutron-star theorists with an 
abundance of astrophysical phenomena to study in more 
detail, and with more precision, than ever before. Among 
these is the formation and evolution of binary pulsars, 
and in particular how this system, and perhaps others 
like it, have formed and evolved to arrive at their current 
configurations. 

We have performed analysis confirming that perhaps 
the channel of binary evolution undergone by this sys- 
tem is somewhat different than that of other double neu- 
tron star (DNS) systems for which orbital parameters and 
component masses have been measured. The remainder 
of this paper will focus on describing the direct observa- 
tional constraints on this evolution that we have gained 
through study of this system. 

FORMATION AND EVOLUTION OF 
THE DOUBLE PULSAR SYSTEM 

It is the evolution of the progenitor to the B pulsar that 
is believed to have caused the differences in current sys- 
tem properties between the PSR J0737-3039A/B sys- 
tem and, for example, the PSR B1913H-16 and PSR 
B1534H-12 DNS binaries as we now observe them [e.g.. 



H. The properties of the PSR B1913H-16 and PSR 
B1534H-12 systems suggest that the companion neutron 
stars (NS) formed via a core-collapse supernova. In par- 
ticular, this is supported by their large respective eccen- 
tricities, e = 0.617 and 0.274, and their high transverse 
velocities, Vtr — 88 and 107kms^' iHUl respectively]. 

The observed parameters of the PSR J0737-3039A/B 
system — e = 0.088 and Vtr = lOkms^' (E] — appear to 
be much lower than typically expected if the second su- 
pernova applied a large natal kick to the system, though 
are not by themselves conclusive evidence against a core- 
collapse supernova event having occuiTed. However, the 
mass of the second-formed NS, PSR J0737-3039B, is 
also low, with = 1.2489±O.OOO7M0. Taken together, 
these properties present important clues that point to a 
different evolutionary path that formed the second NS. It 
has been suggested that this may best be explained by the 
progenitor to pulsar B having gone through an electron- 
capture supernova |3, 7]. This may avoid large supernova 
kicks due to the hypothesized short timescale over which 
this type of event proceeds, which is much shorter than 
that needed for instabilities that produce large kicks to 
develop. Podsiadlowski et al. [3J have computed models 
of this scenario and found that the critical mass for col- 
lapse of an O-Ne-Mg core should range from 1.366 to 

I. 375M0. 

It has been noted by Stairs et al. ItII that the NS 
progenitor must have a minimum mass equal to that 



of the neutron star plus the mass equivalent of the NS 
binding energy. The binding energy of the NS is given 
by ^ O.O84(Mns/M0)2M0 0,%, coiTesponding to 
Mb - 0.13Mq for PSR J0737-3039B, approximately 
equal to the difference between the B pulsar mass and 
the estimates made by Podsiadlowski et al. 1 3] of the pre- 
supernova mass. If indeed this pulsar was formed from a 
collapsing O-Ne-Mg core, this would mean that almost 
no baryonic matter would have been lost in the process. 
The remaining energy will have gone into changing the 
orbital properties of the system, and/or contributing to a 
kick at the time of the supernova. The measured mass 
of PSR J0737-3039B, as well as the low transverse 
velocity of the system, thus present tantalizing evidence 
that the progenitor of the B pulsar may have ended its life 
in this manner 

In addition, it is expected from accretion theory that 
the spin axis of the A pulsar will have become aligned 
with the total angular momentum of the binary system 
(well-approximated by the orbital angular momentum) 
as it gained matter donated by the B pulsar progenitor. 
If the supernova undergone by the B pulsar is close to 
symmetric, this alignment will not be disturbed 1 10]. By 
contrast, if there is a large kick to the system, the result- 
ing misalignment will equal the angle between the orbital 
planes of the system before and after the supernova event 
[e.g.,[il. Several studies [e.g.,[Tl[il[lll3], have exam- 
ined the explosion of the B progenitor Using the timing- 
derived proper motion [6], Stairs et al. |7] predict a post- 
supernova misalignment angle 5 for PSR J0737— 3039A 
of < 11°. Willems et al. [15] predict compatible values 
under different kinematic and progenitor-mass assump- 
tions. If the above studies are correct, measuring a low 
spin-orbit misalignment angle for PSR J0737— 3039A, in 
conjunction with the low system eccentricity and trans- 
verse velocity, would provide crucial evidence in ex- 
plaining the observed properties of this system. 



GEODETIC PRECESSION AND 
LONG-TERM PROFILE CHANGES 

According to general relativity, the spin axis of a pulsar 
in a binary system will precess about the total angular 
momentum vector of the system. This occurs at a rate 
given by lOllil: 
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where in this formulation, nii and are respectively the 
pulsar and companion masses, expressed in solar masses, 
e is the orbital eccentricity, P/, is the orbital period, and 
Tq = GMq/c^ = 4.925490947 ^s is the mass of the Sun 
expressed in units of time. 



The effect of geodetic precession on the spin axis ori- 
entation of this pulsar serves to change our line of sight 
through the pulsar emission region over time. The ex- 
tent to which these effects can be observed depends on 
the spin and orbital geometries of the pulsar system. In 
particular, this is the case for the angle of misalignment 
5 between the pulsar spin and orbital angular momen- 
tum, which defines the opening angle of the cone that is 
swept out by the spin axis over the course of a preces- 
sion period. These geometric parameters can be modeled 
and determined through long-term analysis of the pulse 
profile shape. It has been shown that for PSR B1913H-16, 
the relative amplitudes of the two major pulse compo- 
nents are changing significantly over time, as is the sepa- 
ration between these components flliSl- It has also been 
demonstrated that the PSR B1534H-12 profile displays a 
change in shape over time that is attributed to precession 
by an amount that is consistent with the general relativis- 
tic prediction 11911 . 

PSR J0737-3039A has a precession period of - 75 
years H], much shorter than those of PSR B1913H-16 
and PSR B1534H-12. The effects of geodetic precession 
on the PSR J0737-3039A pulse profile should thus be 
more readily observable on a shorter timescale than for 
PSR B1913H-16 and PSR B1534H-12, provided that 5 is 
non-negligible. 



A NEW SEARCH FOR GEODETIC 
PRECESSION EFFECTS ON PSR 
J0737-3039A 

Observations: past and present 

Using a data set that spanned almost 3 years, Manch- 
ester et al. ||20ll found no evidence that the width or shape 
of the PSR J0737-3039A profile was changing with 
time. They argued that this may be attributed either to 
a small misalignment angle 5 or to the pulsar being at 
or near a precession phase of 0° or 180°. They assumed 
in their analysis that the pulsar emits from a single mag- 
netic pole and so having 5^0° implies a broad beam 
structure that wraps around the neutron star, which seems 
rather unlikely. On the other hand, having the pulsar at a 
special precession phase is statistically improbable. 

To further investigate, we have extended the above 
analysis to include data taken between June 2005 and 
April 2007 with the Green Bank Telescope, using the 
Green Bank Astronomical Signal Processor (GASP) 
backend. These data were obtained through biannual 
concentrated observing campaigns that spanned approx- 
imately 1-2 weeks each, and which typically contained 
5-6 observing sessions that each lasted 6-8 hours. These 
observations were taken exclusively using the 820-MHz 
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FIGURE 1. Pulse profiles and difference residuals over 
nearly two years of GBT observations. High signal-to-noise, 
aligned profiles are shown from each of five concentrated ob- 
serving campaigns. The most recent of these, taken in April 
2007, is shown in the top panel. Each subsequent (shifted and 
scaled) profile and its difference from the April 2007 profile are 
shown and labelled in arbitrary flux units. One can see the lack 
of significant pulse shape change over time. 



receiver, and in general, used up to 16 x 4 MHz frequency 
channels. These data have been flux-calibrated in each 
polarization using the signal from a noise diode source 
that is injected at the receiver. We then calculated an av- 
erage pulse profile for each observing epoch. 

Figure [1] shows the aligned and subtracted profiles for 
each epoch. While there seem to be subtle changes in the 
profile shape over the almost 2 years of observations, it 
is certainly not as apparent as one would expect if the 
spin-orbit misalignment is significantly non-zero. The 
combined Parkes/GBT data set now spans nearly 6 years, 
or - 7.5% of the PSR J0737-3039A precession period; 
the hypothesis that the pulsar existing in a special phase 
of precession thus seems to be an unlikely explanation 
for the observed lack of profile variation. 



Fitting for orbital geometry 

The dependence of the observed pulse width on the 
system geometry can be expressed by the following: 



where 'Pq is half the full pulse width, the angle a is 
that between the spin and magnetic axes, ^ is the angle 
between the A pulsar spin axis and the observer's line of 
sight, and p is the opening angle of the emission cone 
lulll . In using this equation, we are assuming a circular 
emission beam (i.e. no variation in p with latitude on 
the pulsar surface). In their work, Manchester et al. ||20ll 
investigated the effect of using a noncircular beam, and 
found that the results were only marginally affected. We 
thus believe that a circular emission beam is a reasonable 
assumption for this analysis. 

As discussed earlier, geodetic precession causes the 
angle ^, describing our line of sight, to vary over time. 
It can be expressed as a function of time, as well as 
the misalignment angle 5 and epoch of zero precession 
phase Ti, as follows ll22ll : 
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where t is the observing epoch, and / is the orbital in- 
clination. In this analysis we assume that cos / > 0. The 
inclination / = 88.69+0 7g° reported by Kramer et al. 
is sufficiently close to 90 ° that the contribution from the 
first term in equation (O does not change by a significant 
amount by replacing / by 180° — 

If indeed the spin axis of the A pulsar is aligned, or 
nearly aligned, with that of the orbital angular momen- 
tum, then it follows that the spin axis should be nearly 
90 ° from the region of emission. This would mean that 
this pulsar is an orthogonal rotator, presenting the pos- 
sibility that we are viewing emission from both mag- 
netic poles. In investigating the geometry of the PSR 
J0737— 3039A/B system, we performed analysis assum- 
ing both one- and two-pole emission models, where in 
the latter, each pulse component is treated as being emit- 
ted from separate magnetic poles. 



Geometric constraints 

The Parkes and GBT data were obtained at different 
observing frequencies, causing a systematic shift in mea- 
sured pulse width due to pulse evolution with frequency. 
In order to be able to use all the data in a single fit, we 
adjust the measured Parkes widths by subtracting from 
each of those measurements the difference between the 
weighted means of the GBT and Parkes data sets. To ob- 
tain the eiTors on the adjusted Parkes data, the original 
Parkes data uncertainties are added in quadrature to that 
of the difference in the two means. This subtraction is of 
course not physically correct, but it does allow the study 
of pulse shape evolution, since it retains the principle that 
there is zero pulse change within each data set. 

Figure |2] shows the measured pulse widths as- 
suming both one- and two-pole emission from PSR 
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FIGURE 2. Half-profile widths for PSR J0737-3039A as a 
function of MJD for Parkes filterbank data (open circles) and 
GASP data (filled circles). The Parkes data set was adjusted by 
subtracting the difference of the (weighted) average widths of 
the Parkes and GBT data sets. Top: Widths measured at 10% 
of the peak amplitude, assuming single-pole emission. Bottom: 
Widths measured at 25% of the peak amplitude, for the first 
(left) and second (right) pulse components. 



J0737— 3039A. For the former, we take our measure- 
ments at 10% of the peak pulse height for consistency 
with the Manchester et al. [20] analysis. For the latter, 
we re-determined the width of each pulse component 
separately, opting to use the 25% peak flux width for 
these measurements. Below this fractional pulse height, 
there are features in the profile that, when noise is added, 
negatively affect our measurements, making them unre- 
liable. We chose to leave out the initial discovery data 
point from the two-pole analysis, since the extremely 
low signal-to-noise in that profile makes it difficult to 
obtain a reliable width measurement. One can see that 
within each set of data there is no evidence for pulse 
shape change. 

We performed fits on each set of width data using 
equation (|2|i, probing the ;i;^-space over 2 sets of param- 
eter grids. The first pair of values we search are 5 and Ty, 
with < 5 < 180° and with Ti running over one preces- 
sion period centred at MJD 53000. The second grid runs 
overaand5,with0°<a< 180°andO°<5< 180°.In 
all cases, we choose uniform prior distributions for these 
quantities within the respective constraints chosen, as we 
do not have information that would convince us to ex- 
pect otherwise. These fits were done using a Levenberg- 
Marquardt algorithm [e.g., 23], holding each pair of pa- 
rameters fixed at each grid point, allowing the remaining 
two to vary as free parameters, aiTiving at a best-fit 
value correspoding to each {5,Ti) and (a, 5) combina- 
tion, respectively. 

For each grid, we arrived at a joint probability den- 
sity function (PDF), from which we then calculated 
marginalized probabiUty densities for each of the grid pa- 



rameters. The resulting constraints found for a and 5 in 
the one- and two-pole cases are reported in Table [T] In all 
cases we find that a is constrained to be close to 90°, so 
that PSR J0737-3039A is indeed an orthogonal rotator. 
The results shown for 5 are those from the (5, 7\) grid 
fit. We have greater confidence in those upper limits on 
5, for the following reason: There is no reason to prefer a 
specific epoch of zero precession over another. This is a 
quantity that is determined by our line of sight to the pul- 
sar, which we expect a priori to be uniformly distributed. 
In the case of a, however, we know less about the prior 
distribution, and so we choose a uniform prior for lack of 
a more compelling option. 

In addition, the marginalized PDFs for 5 between ° 
and 90° are mirrored by those between 90° and 180°, 
and it has been shown that having 5 < 90° is more 
physically likely, unless the supernova kick is extremely 
large, based on models of the misalignment angle in PSR 
B1913H-16 [24]. We thus quote upper limits based on the 
PDFs for 8 in the region below 90°. We also find T\ to 
be relatively unconstrained, which is expected for a low 
value of 5. It also supports the idea that it has become 
difficult to ascribe this lack of change to the pulsar being 
in a special precession phase, considering that our data 
spans approximately 7.5% of the precession period for 
PSRJ0737-3039A. 

Aberration. Based on our derived geometry for the 
PSR J0737-3039A/B system, we have predicted a 
very small latitudinal abeiTation signal as described by 
Rafikov and Lai [21] (see also [22]), which is expected 
to periodically distort the pulse profile over an orbital pe- 
riod. Using the GBT data, we searched for, but found no 
significant evidence for the effects of aberration on the 
pulse profile of the A pulsar. 



One- or two-pole emission ? 

For each of the one- and two-pole analyses, we cal- 
culated histograms of best-fit values found for p at all 
grid points for the (o:,5) and {8^T\) grid searches we 
performed. In the single-pole emission case, we find 
that it is nearly impossible to avoid solutions with p > 
90°. It seems to be the case that if we constrain PSR 
J0737— 3039A to be emitting from a single magnetic 
pole, that it would have to do so from a fan-shaped beam 
structure extending beyond a single hemisphere of the 
NS, possibly emitted from far in the magnetosphere. This 
would mean that the beam centre is actually on the oppo- 
site side of the pulsar to where we expect, based on the 
profile symmetry. 

In the two-pole analysis, we find that p < 90° is 
clearly favoured, with high occurrences of small values 
of p within that range. In this scenario, it is uncertain as 



TABLE 1. Summary of results for a and S parameter estimation from long-term profile evolution analysis of PSR 
J0737— 3039A. Percentages reflect confidence intervals. 
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to why the pulse peak separation is not 180°. However, 
allowing emission to come from both magnetic poles 
in this pulsar thus avoids the need for exotically large 
beams, or else a reinterpretation of the the location of the 
emission beam centre. 

We thus favour a configuration in which PSR 
J0737— 3039A is an orthogonal rotator with a spin 
axis that is nearly aligned with the orbital angular mo- 
mentum of the system, and that the pulse profile that we 
observe is due to emission from both magnetic poles. 
We then adopt a = 89.6+5 3° ' ™d *e 68.3% upper 
limit 5 < 6.1 °, based on the results from the first pulse 
component, which is the more conservative of the two. 

IMPLICATIONS FOR EVOLUTION 

A low misalignment angle means that it is likely that a 
very small kick was imparted on the system due to the su- 
pernova of the pulsar B progenitor. This agrees with the 
low-kick hypothesis that is favoured by measurements 
finding a low transverse velocity and eccentricity of the 
double pulsar, as well as the low mass of the B pulsar, 
all found from timing measurements of the double pul- 
sar ['6']. It also lends credence to the studies performed 
by Piran and Shaviv 1 13], Willems et al. |14], and Stairs 
et al. |7], which favour a low-mass (< 2Mq) progenitor, 
and a low natal kick (< 100 kms^') given the constraints 
of low space velocity. In the case of the latter study, our 
constraints also agree with their estimate of the misalign- 
ment angle, which they predict to be 0.5 < 5 < 1 1 ° (95% 
confidence). 

This analysis thus supports a scenario in which the 
PSR J0737— 3039B underwent a low mass-loss, rela- 
tively symmetric supernova event. The prominent can- 
didates for such an event are an electron-capture super- 
nova, or a low-mass iron core collapse, and we therefore 
favour one of these scenarios as the one that produced 
the double pulsar system as we now observe it. Discov- 
ery and observations of an increasing number of DNS 
binary systems will help to determine how prevalent this 
type of system is, and will provide further insight into 
this alternate channel of double neutron star formation. 
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